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Abstract 



Hubbard chains with periodically modulated coupling constants in a magnetic field 
exhibit gaps at zero temperature in their magnetic and charge excitations in a variety 
of situations. In addition to fully gapped situations (plateau in the magnetization 
curve and charge gap), we have shown [Phys. Lett. A268, 418 (2000)] that plateaux 
also appear in the presence of massless modes, leading to a plateau with a magne- 
tization m whose value depends continuously on the filling n. Here we detail and 
extend the arguments leading to such doping-dependent magnetization plateaux. 
First we analyze the low-lying excitations using Abelian bosonization. We compute 
the susceptibility and show that due to the constraint of fixed n, it vanishes at low 
temperatures (thus leading to a magnetization plateau) even in the presence of one 
massless mode. Next we study correlation functions and show that one component of 
the superconducting order parameter develops quasi-long-range order on a doping- 
dependent magnetization plateau. We then use perturbation theory in the on-site 
repulsion U to compute the width of these plateaux up to first order in U. Finally, 
we compute groundstate phase diagrams and correlation functions by Lanczos diag- 
onalization of finite clusters, confirming the presence of doping-dependent plateaux 
and their special properties. 

PACS numbers: 71.10.Fd, 71.10.Pm, 75.60.Ej 
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I. INTRODUCTION 



Strongly correlated electron systems in low dimensions are presently a subject of intense 
research. In particular, the magnetism of such systems has revealed very interesting proper- 
ties and it is by now well established that spin chains and spin ladders present plateaux in 
their magnetization curves. It has been shown theoretically that plateaux occur in general 
at rational fractions of the saturation magnetization. The position of these plateaux is sub- 
ject to a quantization condition that involves the volume of a translationally invariant unit 
cell (see e.g. [0-0 )• From the experimental side, different materials have been found which 
exhibit plateaux. One example is a dimerized spin-1 chain || which exhibits a plateau at 
half the saturation magnetization as predicted in ||. A good candidate for a plateau at 
one third of the saturation magnetization is Cu3Cl6(H 2 0)2-2H 8 C4S02 0, though it is not 
yet fully clear whether the proposed frustrated trimer chain model is really appropriate or 
which parameters should be used |jTD[ . 

The most striking examples of plateaux have so far been observed experimentally in the 
materials SrCu 2 (B0 3 )2 [1T| an d NH 4 CuCl 3 JE2] and both constitute again challenges for 
theory. There is generally agreement that SrCu2(BOs)2 is a predominately two-dimensional 
material and how it should be modeled theoretically. Though some progress has been made in 
understanding the origin of some of the observed plateaux ||13|| , it remains a difficult problem 
to compute the complete magnetization process within this model. On the other hand, the 
high-temperature crystal structure of NH 4 CuCl 3 suggests a one- dimensional model, but 
nevertheless it remains unclear what really is the appropriate theoretical model for this 
compound. 

Materials with a ladder structure (see e.g. [14]) are further good candidates for exhibit- 



ing magnetization plateaux. However, since the copper-oxide related materials are strongly 
coupled, plateaux with non-zero magnetization are predicted in a magnetic field range which 
causes difficulties with the present experimental tools. A mechanism yielding plateaux at 
lower values of magnetic fields would therefore be very attractive. As we have shown in 



the case of modulated Hubbard chains |15| , doping may actually provide such a mechanism 
since it allows a continuous variation of the plateau magnetization m with the filling n - 
extending in this particular case also into the low-field region. Doping-dependent magnetiza- 
tion plateaux have recently been also theoretically studied in a different system, namely an 
integrable spin-S* generalization of the t — J chain doped with (S — 1/2) carriers [16|] where, 
however, the appearance of plateaux is restricted to large magnetization values. Another 



example of such a situation occurs in the one-dimensional Kondo lattice model ]T7 where 
unpaired spins behave ferromagnetically thus giving rise to a spontaneous magnetization of 
a value controlled by doping. 

Here we detail and extend our previous study |L5| of the effect of a magnetic field and 
a periodic modulation (p-merization) of the hopping amplitude or the on-site energy on a 
doped one-band Hubbard chain whose Hamiltonian is given by 

L 

H = - ^2t(x) (cl +ha c XjCt + H.c.) + U XXt c *>T c ia c *a 
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h L 



+ ^/i(x) cl a C XyQ - -J2( C l,1°x,1 - C l,l C *d) ■ 



x=l 



Here c x a and c XjCe are electron creation and annihilation operators at site x, a =|, J. the two 
spin orientations and h the external magnetic field. The hopping amplitude t(x) and the 
chemical potential n(x) are taken as periodic in the variable x with period p. 

The one- dimensional Hubbard model with dimerized coupling constants (p = 2) is real- 



ized in a number of real compounds like the organic (super) conductors [18] and the ferroelec 



trie perovskites ]T9|. While some materials in the former class come at quarter filling, one 
frequently finds also realizations of the half-filled Hubbard model. In this case, the model 
( |1.1|) is in the same universality class as a modulated spin-1/2 Heisenberg chain. Realiza- 
tions of the latter exist also at periods p > 2: Some examples of trimerized chains (p = 3) 
have been studied in [|(|^l].[| . 

A technical motivation for resorting to the one-dimensional Hubbard model is that the 
uniform chain is exactly solvable by Bethe Ansatz (BA) for arbitrary values of the on- 
site repulsion U, filling and magnetic field The exact solution can then be used to 
construct a low energy bosonized effective field theory p3|-p5[1 which can in turn be used to 



study perturbations of this model (see e.g. |2B|]). Here we first review some aspects of the 



bosonization description of the Hubbard chain |27| , 28 and extend it for the case of a finite 
magnetic field h ^ 0. 



Focusing on the case of constant chemical potential /i(x) = /i, we have shown in |TJ 



that magnetization plateaux can appear for the model (|1 . 1|) if the density of particles n and 
magnetization m |2£j satisfy 

-(n±m)eZ. (1.2) 

These conditions are commensurability conditions for the up electrons = (n + m)/2 
and down electrons nj_ = (n — m)/2, respectively. More precisely, if both conditions are 
simultaneously satisfied, the system has both charge and spin gaps. On the other hand, if 
only one of these conditions is fulfilled, the filling has to be kept fixed in order to have a 



magnetization plateau. A simple explanation of the conditions ( |L2p can be given |15| in the 
non- interacting limit (U = 0). Then, the Hamiltonian (|1.1| ) can be easily diagonalized and 
is found to have p bands e x (k) (see section [TT1| for more details). The magnetic field breaks 
the symmetry between up- and down-spin electrons by shifting their chemical potentials by 
opposite amounts. It is then possible that one chemical potential (say for the up electrons) 
lies in one of the p — 1 band gaps while the other (for the down electrons) is in the middle of 
a band. This situation leads to a doping-dependent plateau, if one imposes the constraint 
of fixed filling n (and only in this case). Then the magnetization can be increased only by 
moving an electron from the down-spin band into the up-spin band which requires a finite 
energy or equivalently a finite change of magnetic field, leading to a plateau. However, since 
the filling of the down-spin electrons remains adjustable, one obtains a doping-dependent 
value of the magnetization at the plateau. 

Finally, we have also shown [Tj| that a charge gap opens if the combination pn 6 Z of 
the two conditions ( |1.2| ) is satisfied. The latter case generalizes the well known charge gap 
at half filling (n = 1) as well as the charge gap at quarter filling in the dimerized Hubbard 



chain (n = 1/2, p = 2) |30|-f£|. 
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The plan of this paper is as follows: In subsection |II A| we briefly review the bosonization 
approach for the Hubbard chain for arbitrary filling and on-site Coulomb repulsion U in 
the presence of an external magnetic field (our conventions are summarized in appendix [A] 
and details on the bosonization approach in appendix ||) |23|,|25| . In subsection |IIB| we then 
use this bosonization scheme to study the effect of a modulation of the hopping amplitudes 
and the on-site energy n(x) and find the conditions under which a plateau is present. The 
appearance of plateaux for irrational values of the magnetization and superconducting cor- 
relations are analyzed in subsections U. C| and II D| , respectively. In section |T| we study 
the limit of small U perturbatively and show that the doping-dependent plateaux are also 
present there. Then we study the groundstate phase diagram (subsection [IV A] ) and corre- 
lation functions (subsection [IV B| ) numerically on finite size systems by means of Lanczos 
diagonalization. Finally, we summarize our results in section M, discuss some experimental 
settings where the features presented in this paper could be observed and point out open 
routes for further research. 



II. BOSONIZATION APPROACH 

A. Field theory description of the Hubbard chain in a magnetic field 

In this section, we summarize the analysis of the Hubbard model in a magnetic field 
using Abelian bosonization. For further details see Ref. [25|. The lattice Hamiltonian is the 
standard one, i.e. (|1.1|) with constant t(x) = t, /i(x) = /i: 



(2.1) 



This model has been exactly solved by BA already in 1968 |22fl but it took until 1990 for the 
correlation functions to be computed by combining BA results with Conformal Field Theory 
(CFT) techniques ||23|| . Spin-charge separation is a well known feature of the Hubbard chain 
at zero magnetic field. Interestingly, it is no longer spin and charge degrees of freedom that 
are separated if an external magnetic field is switched on p3 |. Nevertheless it has been 
shown that in the presence of a magnetic field, the spectrum of low energy excitations can 
be described by a semi-direct product of two CFT's with central charges c = 1 [23]. This 



in turn implies that the model is still in the universality class of the Tomonaga-Luttinger 
(TL) liquid and therefore allows for a bosonization treatment. 
In order to proceed, we write the fermion operator as 



~ e 



1pR >a ( X ) +••• 



(2.2) 
(2.3) 



where &f )Q are the Fermi momenta for up and down spin electrons and are the chiral 

components of two bosonic fields, introduced as usual in order to bosonize the spin up and 
down chiral fermion operators ipR,L,a- (Our conventions are settled in appendix |A]). The 



4 



dots stand for higher order terms some of which are written explicitly in appendix [B| They 
take into account the corrections arising from the curvature of the dispersion relation due to 
the Coulomb interaction. For non-zero Hubbard repulsion U and magnetic field h, the low 
energy effective Hamiltonian corresponding to ( [2.1|) written in terms of the bosonic fields 0j 
and 4>i has a complicated form, mixing up and down degrees of freedom. The crucial step to 
obtain a simpler bosonized Hamiltonian is to consider the Hamiltonian of a generalized (two 
component) TL model and identify the excitations of the latter with the exact BA ones for 
the model ( |2.1| ), providing in this way a non-perturbative bosonic representation of the low 
energy sector of the full Hamiltonian ( |2.1| ). This program has been carried out in Ref. 
and we just quote here the final result. The fixed point (i.e. neglecting all irrelevant terms) 
bosonized Hamiltonian reads 



H 



dx 



u. 



u. 



(2.4) 



where = (4>r,-\, 4>l,-\, 4>r,i, <A&,j.)- The matrices A c , s have the following form: 



Ac 



/ a c , s + b, 



c,s "c.s 



c r c d r 



J C,S '"C.S 



'C,S 



a c .s + b c 



+ d c 



Cc.s d c 

- d c 

+ fc 
-fc 



dr S \ 



+ d 

+ U. j 



c,s 
c.s 



(2.5) 



where 



d s 



-1\2 



Z sc (Z C c — Zsc) 



1\2 



a s — (Z sc 



b c — (Z cc ~ Z. ^ 
e c = (Z- 1 + Z, 



sc I C c 
cs I 



[z- 1 + z- 



Zc^iZcc + Z cs l , 



fc 



zi 



2 



Z S c(Z ss 



f 



ss 



(2.6) 
(2.7) 



In these expressions Z^ (resp. Z^ 1 ), i,j = c, s, are the entries of the dressed charge matrix 
Z (resp. its inverse Z~ l ) taken at the Fermi points 



Zcc Z cs 
7 7 

^ sc ^ ss 



(2.8) 



These matrix elements are solutions of a set of coupled integral equations obtained from the 
BA |23|] and depend on the coupling U, the chemical potential \i and the magnetic field h. 
They can be in turn related to physical thermodynamic quantities [^3 . 
Substituting for the bosonic fields 



7 7 

^ SS " ss 

det Z V Z sc Z sc 



(2.9) 



and for their dual fields 



^ sc *-* sc 

ZcS Z BC 



0, 



(2.10) 



the Hamiltonian takes the form 
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J2^[dx [(d^f + [dM] , (2.H) 

i=c,s 

where (ft = 4>r + 4>l and 6 = (pR — 0^. 

At zero magnetic field, the matrix Z reduces to 

Z <* = °>=( { /2 1/W ■ < 2 - 12 > 

with £ = [/). In this case we recover the well known expressions for the charge and spin 
fields 

<Pc = ~ (0 T + 0l) , 0* = 4= (0T - 0l) > (2-13) 

where the compactification radius of the spin field {i.e. the parameter which indicates the 
period of <p s , <fi s = 4> s + 2nR s , R s = l/y/2n) ]34| is fixed by the SU(2) symmetry of the 
spin sector. The radius for the charge field, on the other hand, depends on the chemical 
potential /i and the Coulomb coupling U. Furthermore, for h = the charge and spin 
degrees of freedom are completely decoupled. 

It should be noted that for m ^ 0, the fields arising in the diagonalized form of the 
bosonic Hamiltonian (|2.11|) are no longer the charge and spin fields even though they have 
been labeled by 'c' and L s\ For example, the charge field is in general given by <fi^ + 0j_ = 

Zcc4>c — Z cs 4> s . 

For generic values of the parameters of the model ( j2.1|) , we can now write down for 
example the bosonized expression for the charge density operator: 

p{x) = ip\ipi(x) +i/;lip l (x) 

= —i=d x (Z cc (j) c - Z cs (j) s ) + a p sinf/c+x - ^/n (Z cc (/) c - Z cs (/) s )] 
V 71 " 

x cos[/c_x - ^/n ((Zee - 2Z SC )0 C - (Z cs - 2Z ss )(p s )] 
+ 6 p sin(2A; + a; - V4^(Z cc <f) c - Z cs <f) s )) , (2.14) 

where a p , b p are non-universal constants, whose numerical values are known only in special 
cases. Details on how such expressions are obtained are given in appendix [B]. Formulae of 
the type (|2.14 ) are our fundamental bosonization rules. 



B. Space dependent modulations 

In the present subsection we study two different perturbations of the Hubbard chain 
( p.lD , which consist in space dependent modulations of certain parameters. In particular 
we shall consider a space dependent modulation of the hopping amplitude t(x) and of the 
on-site energy [x(x). 
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1. Modulated hopping amplitude 



In this case, the Hamiltonian reads as in eq. (IT]) with /i(x) = const and t(x) = t if 



x 7^ Ip and t{lp) — t' — t + 5, with p, I integers and p fixed. This is equivalent to the uniform 
Hubbard Hamiltonian (|2.1|) perturbed by the term 



H pert = -5 ^ (Cy.qCV+l.Q + H.c.) . (2.15) 

x'=lp,a 

At half filling and for large U, a standard second order perturbative computation in 1/U 
shows that the effective Hamiltonian is given by 

# = £^4-4+1, (2.16) 



X 



thus leading to the p-merized Heisenberg chain studied in 0. It was predicted there that 
magnetization plateaux occur when the condition |(1 — m) G Z is satisfied p9| . We use now 



Abelian bosonization techniques to analyze the more general case of the model ( |1 . 1| ) in the 
small 5 (weak p-merization) limit. 

Using the bosonization dictionary given in appendix [B], we find the expression for the 
continuum limit of the lattice perturbation (|2.15|) 



Open = Ai sin[A; + /2 + pk + x - \pH (Z cc (f) c - Z cs (f) s )} 

x cos[A;_/2 + pk_x - ^ ' {{Z cc - 2Z sc )<f) c - {Z cs - 2Z ss )<f) s )] 

+A 2 sin[fc + + 2pk + x - V^(Z cc (f) c - Z cs (p s )] , (2.17) 



where Ai, A 2 oc S and the Fermi momenta are k + = k F ^ + k F ^ = nn, fc_ = kp^ — kp,i = ^m, 
where n is the filling and m is the magnetization. The presence of a factor p in the oscillating 
part will play an important role in the following. 
The operator 

A 3 cos[A;_ + 2pk^x - 2^((Z CC - 2Z sc )<p c - (Z cs - 2Z ss )(j) s )} , (2.18) 



with A3 oc 5 2 is radiatively generated from the first term in (|2.17 ) and must therefore be 
included as well. 

In the case of zero magnetic field the dressed charge matrix is given by ( |2.12j ) and we 
have then a neat separation between charge and spin fields. The most relevant perturbation 
takes the form 

Opert = Ai si 11 ! - ^ - + pnirx — \/7r£0c] cos[V2iT(j) s ] + A2 sin[7m + 2pmrx — VAtt^c] . (2.19) 

The marginal operator associated with A3 contains only the spin field, its dimension (fixed 
by the SU(2) symmetry) is 2 and it is marginally irrelevant. A term like this is already 
present in the original model and is also marginally irrelevant. Hence, for S small enough, 
this term can be absorbed in the original marginally irrelevant perturbation term without 
changing its relevance character. 
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The A2 term affects only the charge degrees of freedom and its dimension runs from 1, 
for U — > 00 to 2, for U = 0, being then always relevant for the cases of interest. We can 
therefore conclude that the charge field is massive whenever this operator is commensurate, 
which in turn happens if the condition pn G Z is satisfied. 

If this happens, we can integrate out the massive charge degrees of freedom which leaves 
us with an effective theory for the spin degrees of freedom. This effective theory is massless 
except when the operator associated with Ai becomes also commensurate, i.e. if the condition 
^ G Z is satisfied. In that case, we have also a spin gap in the system. 

These considerations are easily generalized to the case of non-zero magnetization as 
long as the condition pn G Z is satisfied. In this case, the A2 term in ( [2.17 ) is always 



commensurate. Since it contains only the proper charge field fa + a charge gap opens 
for all values of m at these commensurate values of the filling. The condition pn G Z is also 
satisfied when the two conditions ( p..2|) are simultaneously satisfied. In this case, also the Ai 
term in ( 2.19Q becomes commensurate thus leading also to a spin gap. 



In particular for p — 2, 3, we predict the following fully gapped situations: 
p = 2: Half filling (n = 1): gap for the charge, and plateau for m = 0. 



Quarter filling (n = 1/2) (and also n = 3/2): gap for the charge [pOH 33||, and plateau 
for m = ±1/2. 

p = 3: n = 1, n = 1/3 and n = 5/3: gap for the charge, and plateau for m = ±1/3. 
n = 2/3 and n = 4/3: gap for the charge, and plateau for m = ±2/3, 0. 



The final case where only one of the conditions ( |1.2| ) holds is more complicated since then 
the charge and spin degrees of freedom can no longer be separated. We therefore postpone 
discussion of this case a little. 



2. Modulated on-site energy 

Now we consider the Hubbard chain ( |1 . 1|) with a uniform hopping amplitude t(x) = t 
but a periodic modulation of the chemical potential fi(x) = fi if x 7^ Ip and fi(lp) = \i + 5fi, 
with p, I integers, p fixed. This is equivalent to the uniform chain Q2.1|) plus an on-site energy 
term which reads 

x'=lp,a 



The case p — 2, h — has been studied in detail in p5| . 

In the continuum limit the perturbing operator ( 2.20Q becomes: 

0' per t = Ai sm[pk + x - y/n (Zccfa - Z cs (f) s )] cos[p/c_x - ^tx((Z cc - 2Z SC )4> C - (Z cs - 2Z SS )4> S )] 
±A 2 sm(2pk + x - V4n(Z cc (j) c - Z cs (p s )) , (2.21) 

which radiatively generates a term of the form 

A 3 cos[2pA;_x - 2^{{Z CC - 2Z sc )fa - {Z cs - 2Z ss )fa)] . (2.22) 



S 



The only difference with respect to the previous case, eqs. (|2.17|) , (|2.18|) , is a phase in each 
of the sines or cosines, which however plays a role only at half filling and zero magnetic field 
(see also J35|). Apart from this particular case, the conclusions remain the same as in the 
previous case. 



C. Partial gap: irrational plateaux 

We have shown in the previous subsection that, when both commensurability conditions 
( |1.2|) are satisfied, the spectrum of the model (|1.1|) is fully gapped. It remains to understand 
what happens if only one of these conditions is satisfied. In this case, apparently one degree 
of freedom remains massless. We will show that the system nevertheless exhibits a gap for 
magnetic excitations (i.e. excitations changing the value of the magnetization), provided 
the total charge (i.e. the filling n) is kept fixed. 

For the sake of simplicity we will restrict ourselves to the dimerized chain (p = 2) in this 
subsection though the argument can be generalized easily to p > 2. Suppose that 

p 

— (n — m) = n — m = 2n± E Z , (2.23) 

but that n + m is not an integer, i. e. the commensurability condition is fulfilled only for 
down electrons. Assume first that there is no interaction between up and down electrons 
(U = 0). Then we are led to analyze the excitation spectrum of the following Hamiltonian 
in a system of length L: 

H = Jdx^ [(d x ^f + (OA) 2 } + | [(d x <Pif + (dA) 2 } + A oos(2v^) , (2.24) 
with total magnetization 

M = [ dx -^d x (^ - = -^(0 T - ^)\o ■ (2-25) 
Jo V 71 " V 71 " 

Motivated by experimental realizations of Hubbard systems, where typically doping is fixed, 
we also impose the constraint that the total particle number is fixed: 

N = I dx ^=d x (^ + 4>i) = ~^(h + <f>l)\o ■ (2-26) 

'TV \ TV 



From eqs. ( |2.25| ) and ( |2.26| ) we see that the fields <p^i satisfy the following boundary condi- 
tions: 

24>^ = ^(N + M) , (2.27) 

20|lo = V^(N-M) . (2.28) 

Notice, furthermore, that the fields are compactified, i.e. they satisfy the periodicity 
condition 

<h,l - h,l + v 7 ^ • (2-29) 
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Therefore, in a semiclassical picture, the vacuum configuration for 0j is 



M x ) = irr( N + M ) x + const ■ ( 2 - 30 ) 
On the other hand, for <f>^ the vacuum configuration is a kink 

</> l (x)=k(x) , (2.31) 

where k(x) is a configuration interpolating between two minima of the cosine potential in the 
Hamiltonian ( |2.24| ) and satisfying the boundary condition ( 2.28 ). Now we change the total 



magnetization keeping the total number of particles fixed. The lowest energy excitation 
of this type consists in reversing the spin of a particle, which corresponds to the change 
M — > M + 2. The new boundary conditions become then 

20 T |£ = ^(N + M + 2) , (2.32) 
20ilo = V*(N-M-2) . (2.33) 

The new vacuum configuration for </h is therefore 



T ( x ) = (N + M + 2)x + const , (2.34) 

and it is straightforward to show that the difference in energy with respect to the original 
one is linear in 1/L. On the contrary, changing the configuration of the kink requires a 
finite amount of energy (proportional to A) because the new configuration is in a different 
topological sector. This corresponds to the presence of a gap in the spectrum of magnetic 
excitations, and therefore of a plateau in the magnetization curve. 

To support the previous conclusion further, we analyze the magnetic susceptibility for 
a chain of finite size L. It is given by the integral of the correlation function of the spin 
density operator -^d x (4>i — 4>i). Since the down sector is gapped it does not contribute to 
the zero temperature limit of the susceptibility. Let us therefore focus on the up sector, 
which is apparently massless but constrained to be in a particular topological sector. One 
can easily see that determination of the susceptibility amounts to calculate 

dx d x A \ . (2.35) 




o 



For the free massless sector, the Hamiltonian in a finite size L can be written in Fourier 
space for each topological sector as (see for example 



1 



q 2 (f>- q (p q + K q 2 e_ q e q 



where Q and J stand for the particle number and current zero modes: 
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and the summation over q is for oscillatory modes. If the global constraint is not present, 
one has to sum over all possible values of Q, i.e. one has to compute 



X 



i f^Tr (eM~PH L )Q 2 ) - (lie (exp(-pH L )Q)j j . (2.37) 



The local part (the oscillator modes) decouples, and if the constraint is not imposed we 
obtain the standard result 




X 



If we now impose the global constraint on N, due to the gap in the down sector as discussed 
above, all the sectors will be exponentially suppressed, except the sector Q = (JV + M)/2. 
Therefore, for small enough temperature, the distribution exp(— P^b-) has to be replaced 
by a delta function in Q = (N + M)/2, giving 

X = ^ Z ((Q 2 )-(Q) 2 )=0. (2.39) 

We find then an exotic situation in which we have simultaneously algebraic decay of correla- 
tion functions, since the local dynamics is massless, but zero magnetic susceptibility, due to 
the global constraint imposed on the system. The only somewhat similar situations we are 
aware of include plateau states of strongly frustrated spin ladders with gapless non-magnetic 
excitations [§,37,38] as well as the large number of singlets inside the gap of the Heisenberg 



antiferromagnet on a Kagome lattice . 



D. Superconducting fluctuations 

Having found a situation with a gap which can be attributed to magnetic excitations and 
another massless degree of freedom, one may wonder whether superconducting fluctuations 
develop. We therefore now briefly analyze the correlators of the superconducting order 
parameter. In the presence of a magnetic field, the superconducting order parameter has 
four components which read on the lattice: 

Aq, )/3 = c x+1)0t c X)l3 . (2.40) 

For h = 0, these components can be grouped in a triplet (t) and a singlet (s). On the lattice, 
the corresponding S z = components can be chosen as 

Kf = c*,T c *+U ± c x ,iC x+1A . (2.41) 
In the continuum, using eq. ( |Alj) this leads to the following expression 

A M = e- ifc -^ R>T ^ La (e ifc ^ =F e"^) + e^L^i (^ kFi T <^ FT ) 

+ e~ ik+x ^R,^R A (e~^ =F e"^T) + e lk + x ^ L ^ LA (e lk ^ t e lk ^) . (2.42) 
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In particular, for zero magnetic field k F "\ = kpi = k F and neglecting "2fcp" terms, eqs. (|2.42|) 
reduce to the standard ones (see for example Ref. 2q|): 



At = 2i sin k F (i) R ^il) L>l + ipR,iip L ,-\) , A s = 2 cos k F {^R,^L,i - i^r^l^x) . (2.43) 
For general m, one finds instead 

A t ~ -L e V*(*T+»i) cos( v / ^(0T - 0i ) - , A * ~ -e l ^ (flf+9l) sin(v^F(0 T - </>,) - . 
na 7ra 

(2.44) 

These expressions show that the correlators associated to the order parameters A t and A s 
decay exponentially, even in the partially massless plateau phases. Indeed, since the S z = 
components of A a)j g are products of up and down degrees of freedom, it is sufficient that one 
of them is gapped in order to lead to an exponential decay of the composite object. 
On the other hand, the diagonal components are bosonized as 

A a , a ~ 2 cos(k F , a )e iV ^ ea + e -ik F , a (i+2x) e iV^(<t> a +e a ) + e ik F>a (i+2x) e iV^(-4> a +9 a ) + _ ? ^A5) 

where the dots include terms which mix <f>^ with <p^. It is then clear that on a doping- 
dependent plateau, where only one of the fields is gapful, only one of the correlators 
(A^ aa A a>a ) decays exponentially, but the other exhibits algebraic behavior. In fact, all 
fields involving only the gapless spin component decay algebraically. In particular, also the 
two-point correlator of c x>a decays algebraically if A aya exhibits quasi-long-range order. The 
algebraic decay of the latter should therefore not be taken as a sign for superconductivity, 
but is interesting nevertheless. 



III. SMALL- [/ LIMIT 

The previous section was dedicated to the bosonization approach to the p-merized Hub- 
bard chain in the small p-merization limit. In the present section, we give a further argu- 
ment for doping dependent plateaux, valid in the low U limit but at arbitrary p-merization 
strength. For the sake of simplicity, we will concentrate on the case of modulated hopping 
amplitude t(x) = t' for x a multiple of p, otherwise t(x) = t and constant p,(x) = fi, but the 
arguments can be generalized easily. 

First we diagonalize the Hamiltonian ( |1.1[ ) at U = by a unitary transformation 

L /p p 

dk,a = ~7f ^ a k,j C xp+j,o- ■ (3-1) 

x=l j=l 

In order for the kinetic part of the Hamiltonian (|1 . 1| ) to take the form 

A=l a 

the coefficients have to satisfy the following eigenvalue equation: 
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/ t ••• t'e- ik \ 

tot'-. 
o t •■. '•. '•. ; 
; ■•. ■•. ■•. t o 
o '-.tot 

\t'e ik •■■ t / 

The resulting p energy bands e x (k) are illustrated in Fig. 2 of for p = 2 and in Fig. 2 
of [0 for p = 3 (note that in the latter case, the energy e was plotted with the wrong sign 
which can be absorbed by shifting k — > k + 7r) . 

In the sequel we first work out the simpler case p = 2 and then generalize to p > 3. 



l fc,2 



(3.3) 



A. Case p = 2 

In the dimerized case, the eigenvalue equation (CO) reduces to 



The eigenvalue problem (|3T4]) is solved readily, yielding 

e^Jfc) = ±Vt 2 + t' 2 + 2tt> cosfc, 

a f = z F t/t±*f!!L, (3.5) 
1 V t + t'e lk v ; 



t + £'e 
£ + t'e- ife ' 

The inverse of the transformation (B.ll) is: 



fc 

Eigenstates |{^ CT }) of the free Hamiltonian Hq are now written down by simply specifying 
the momenta k* a occupied in the various bands. Now we treat the Coulomb repulsion 

L 

Hi = U^2n xA n xd (3.7) 

x=l 

in first order perturbation theory. 

To proceed further, we assume that none of the bands are half-filled. Then Hj has only 
diagonal terms (i.e. Umklapp scattering is absent) which are readily evaluated as (denoting 
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u 



({kl}\Hj\{kl}} = y({^JlEKT + <T)K"M + <i)l{^}) 



L 

U 

T 

u 

1 



k,k' 



(N~ + N+)(N~ + N+) 
L(n^ + nf)(nT + nt) . 



(3.8) 



In the second line, we have defined the number of particles with spin a in band A by 
= pUf . The densities have been normalized such that — 1 for a completely 

filled band. 

Similarly, expectation values of the number operators Y^=\ n x,a gi ye L (n~ +n+) /2. 
Putting everything together, we find the energy of the Hamiltonian ( |1 . 1| ) to first order in U 

as 

U 



A k x 



(3.9) 



Assume now that the '— ' bands are both partially filled. Then Q3.9| ) specializes to 

n CT 7r 

E / L = -^Y. J ^W + jn t »l + fh+»l)-^t-«l)» (3-10) 



where n a = n u Setting n = + and fixing nj ~ 1, we find from the condition that 
it does not require energy to flip | to j, spins, 



h Cl = e 



(n)-e-((n-l)n) + u(^-l 



(3.11) 



On the other hand, if we consider the case of a completely filled "— , j" band and partially 
filled "+, |" and "— , j" bands, eq. ( |3.9|) specializes as follows: 

E/L = i- J dke + (k) + ±- Jdke-(k) + ±- J dke-(k) + j(l + n+) ni 



(l-n+)7r U 

+ |(1 + n+ + n ; ) - |(1 + n+ - raj . 
Setting n = 1 + + and fixing ~ 0, we find in the same way as before that 



h 



,,, - e + (7r)- e -((n-l)7r) + [/ Vo 



1 



(3.12) 



(3.13) 



Using eqs. ( 3.1 3|) and ( 3.11 ) we find that the width of the plateau at fixed N is not affected 
by the on-site Coulomb repulsion to first order in U: 
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h C2 - h Cl = 2 \t - t'\ + 0{U 2 ) . (3.14) 

This ensures the presence of a doping-dependent plateau with m = 1 — n in the low U limit. 
The absence of a first-order correction in U to the width in (|3.14 ) can be traced to the 



mean-field form Q3.8|) of the matrix elements of the on-site repulsion Hi. This in turn is due 
to the fact that \af | = 1 for all k as can be seen from ( |3.5| ) and is a manifestation of the 
symmetry between the lower and upper band. Both the mean-field form of the interaction 
as well as the absence of a first-order correction to the plateau width are special properties 
of the case p = 2, as will become clear in the following discussion of the case p > 3. 



B. Case p > 3 



For general p, the diagonalization ( |3.3| ) is more complicated leading to the absence of ex- 
plicit expressions such as (|3.6j ). Nevertheless, we can still use unitarity of the transformation 
( p\ip to formally invert it 



^EE e ~ 



ikx *X j\ 
a k,j a k,cr ■ 



First, we look at the transformation of number operators n x ^ = c xa c x ^ a 



n 



fc,<7 



(3.15) 



i p 

E n ^ = -EE 



n. 



k,X j=l 



E n t- 

fc,A 



(3.16) 



Here, we have just used that Y^j=i \ a k,j\ 2 = P- 

The diagonal terms of the interaction (3/7) can now be treated similarly as for p = 2. 
Instead of (|3~8|) one finds for general p 



y 3=1 k,\ 
11 P 

=—Y y 

F 3=1 



k,j I 



A I 2 
l k ,j I 



'.T E \ a k',j U k' ,l\{kj,a}) 



fc',A' 



E 



\n x ' 
\ a k',j 



2>T 



fc,A 

occupied 



(3.17) 



Next we pass to the thermodynamic limit which leads to replacing sums by integrals. Due 
to (|3.16| ), integrals and differentials over densities can be replaced by integrals in fc-space. 
We work at a fixed particle number n, which implies 



= dn = dn^ + dn 



dn\ 



-dn^ . 



(3.18) 



Now we concentrate on the situation where all bands A < Ao are completely occupied with 
up spins and those with A > Ao do not contain any up spins generalizing thus the reasoning 
of the previous section. The band A' is partially occupied with down spins, those with 
A' < A' are completely filled with down spins while those with A' > A' do not contain any 
down spins. For a partially filled band v denote finally the occupied states by \k\, k£\. 
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Then we can generalize ( |3.10|) to first order in U as follows [41 



E/L = u + v + ^ (n T +nj.) - ^ (n T - n{) 



with 



u 



2np 



A<A, 



/ dke x (k)+ / dk'e x '°(k')+ J2 / dk'e x \k') 



and 



47T 2 p 3 



EE 

j=l a<A 



dk \ a 



A |2 
k,j\ 



fc A ° 



This yields for the lower boundary of the associated plateau 



h Cl = e 



A„ 



{k^)-e\k^)+U{\ ) + O(U 2 ) 



with 



U{v) 



U 



2vrp 2 



J2\ a k,j\ 2 

v 

-EE 

j=l A<A 



£- A ° 



fc A ° 



dk' 



E 



eifc' 



A |2 



u „• 



For the corresponding upper boundary one finds 

fr C2 = e A o+ 1 (^' +1 ) - e x '°{k$) +U{X + l) + (U 2 ) 
Eqs. dp|) and (ggp imply that 



^C2 ^Cl f 



with 



27rp' 



E 

i=i 



7 A +1 
Vo +1 



l X ° 



fc A ° 



dk' 



x 'o 



A'<A' 



dk' 



a k',j 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



The first-order contribution does in general not vanish (for p > 2). It can be estimated as 
1 



*fc A ° 



A'<A' 



P 



(3.27) 



which shows that in principle it can be of order one. 
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FIG. 1. Value of the first-order correction Oi to the plateau width as given by ( 3.26p for p = 3 
and n-f = 2. 



In general, it is not difficult to evaluate the first-order contribution ( |3.26|) to the plateau 
width numerically. We will illustrate this now for p = 3. First notice that, for the case 
p = 3, the contribution from the kinetic energy can be readily evaluated as 



V8t 2 + t' 2 - 3f 



Now we fix 7ij : 
one has that kf, 



2, i.e. 

= k 2 = 



the lowest two bands of up electrons are completely filled 
in ( EQ6 



(3.28) 
Then 

and evaluation of 



Numerical diagonalization of (|3 
the remaining integrals in ( |3.26|) then leads to Fig. |l[ Note that in the conventions of the 
other sections (where < n < 2) this corresponds to the plateau with m = 4/3 — n. The 
numerical data satisfies 0\ — > —0\ as nj_ — > 3 — ri|. This implies in particular that the values 
of 0\ can be both positive and negative, corresponding to an enhancement or reduction of 
the plateau width, respectively. Furthermore, for t — > and rij < 1, the linear behaviour of 
( P-27| ) is reproduced, though with a coefficient which is 1/6, i.e. by a factor of two smaller 
than in the estimate. The maximal values attained are ±1/6 for nj, — > 1 or 2, respectively 
and t — > 0. This shows that the doping- dependent plateaux should be stable features for 
p = 3 as well. 

We conclude this section by noting that the calculations are also valid for the on-site p- 
merized energy. The free Hamiltonian Ho to be diagonalized is modified, but the conclusions 
remain qualitatively unchanged. 
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IV. LANCZOS DIAGONALIZATION 



Finally, we have performed Lanczos diagonalizations of the Hamiltonian ( |1 . 1[ ) for p = 2 
and p = 3 with constant p(x) = p and periodic boundary conditions on finite lattices 
in order to further support the previous results. The particle numbers nj and nj. have 
been used as quantum numbers and translational symmetry was exploited. Furthermore, 
reflection symmetry was exploited for k = 0,7r and spin inversion for — nj_. 



A. Groundstate phase diagrams 

Computations have been performed mainly for one choice of parameters due to the large 
number of sectors for which the groundstate energy had to be found (for p = 2 and L = 16 
of the order of 10 3 sectors). Keeping (|2.16|) in mind, we have chosen the parameters U = 3t 



and t' = 0.7t in order to look at a situation sufficiently different from the limiting cases 
discussed before, i.e. both intermediate (to strong, as compared to the band-width) on-site 
repulsion U and intermediate t'/t. 

For the interpretation of our Lanczos results to be presented below, it is useful to re- 
member the following consequences of particle-hole symmetry on a finite size lattice (see 
2"2" , 32 ] and references therein): For L even, the groundstate phase-diagram of the Hub- 



bard chain with periodic boundary conditions is symmetric under p i— ► — U — p (with our 
conventions), while for L odd the particle-hole transformation interchanges periodic and 
antiperiodic boundary conditions. 

Our numerical results for the groundstate phase diagram are shown in Figs. ||] and § for 
p = 2 and p = 3, respectively. The polygons in the figures denote regions in the (/x,/i)-plane 
where the groundstate has a fixed filling n and magnetization m at the given system size L 
(those values of m and n which are common to all investigated system sizes are indicated in 
the figures). The groundstate phase diagrams are symmetric under spin inversion (m — > — m 
when h — > —h) and as mentioned before, for even L also under a particle-hole transformation 
(n — > 2 — n when p — > —U — p). Therefore, for even L we show only the quadrant with 
H > —U /2, h > and for odd L only the region with h > 0. 



The schematic groundstate phase diagrams in |15j were in fact based on parts of these 
results and the reader may wish to use them as a guide to the diagrams at finite size. 

We note that for the saturated case n± = (and by particle-hole symmetry also for 
= 1), the Coulomb repulsion is not effective and the non-interacting result {e ± (k) given 
by eq. ( |3.5| ) for p = 2) can be used to determine the transitions between different particle 
numbers. Complete agreement between this analytical computation and the corresponding 
numerical results in Figs. and [5] is found. This also guides the interpretation of the finite- 
size data since it follows in particular that the completely gapped situations at saturation 
are those with pn 6 Z in the thermodynamic limit. Such a guide is useful since the fermions 
behave differently for even and odd particle numbers thus yielding non-monotonic finite-size 
effects which can be still strong for the small systems sizes considered here. In the particular 
case of the gapped states with pn G Z at m = n (or m = 2 — n), the corresponding 
groundstate always has an even number of fermions when L/p is even while cases with 
an odd number occur when L/p is odd. This leads to vanishing finite-size effects for the 
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(I/t 

FIG. 2. Groundstate phase diagram of the dimerized chain (p = 2) with U = 3t, t' = 0.7t. 
In a) the lines are for L = 6 (dotted), L = 10 (dashed) and L = 14 (full) while in b) they are for 
L = 8 (dotted), L = 12 (dashed) and L = 16 (full). 
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FIG. 3. Groundstate phase diagram of the trimerized chain (p = 3) with U = 3t, t' = 0.7t. 
In a) the lines are for L = 9 (dashed) and L = 15 (full) while in b) they are for L = 6 (dotted), 
L = 12 (dashed) and L = 18 (full). Note that the L = 18 data in b) is incomplete for n > 2/3 
(e.g. only m > 1/3 for n = 1). 
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transition lines in the former case, but not always in the latter. Since even and odd L/p 
behave differently, we show separate figures for the two cases. 

For p = 2, one can quite clearly recognize the fully gapped situations at (n, m) = (1,0), 
(1, 1), (1/2, 1/2) and (0,0) from the finite-L data shown in Fig. |2|. Also the charge gap at 
half filling (n = 1) is obvious. The most interesting region is the doping- dependent plateau 
with m = 1 — n which is a stable feature in Fig. ^|a), but less clear in Fig. 0b). Still, in 
the latter case the region of stability of states with m = 1 — n can be seen to increase with 
increasing system size, thus supporting the presence of a gap. Just the charge gap at quarter 
filling (n = 1/2) is not distinct in this numerical data, but it is known to be small anyway 
for these parameters [B(J. 

The case p = 3 is shown in Fig. ||. There is clear evidence for the expected fully gapped 
situations (the labeled regions in the figure) as well as the charge gap at half filling. There 
is also evidence for the charge gap at n = 2/3 and the equivalent case n = 4/3, just the 
charge gap at n = 1/3,5/3 is again difficult to see. Also the expected doping-dependent 
magnetization plateau with m = \n — 2/3 1 can be recognized in Fig. |3]a). By particle- hole 
symmetry, the plateau with m = |4/3 — n\ must be present as well though it is more difficult 
to recognize. The finite-size behavior of its stability region in Fig. |3]a) (and of both plateaux 
in the case of Fig. |3|b)) can again be taken as an indication that it will indeed be present in 
the thermodynamic limit. 



B. Correlation functions 

Having provided also numerical evidence for the existence of doping-dependent magne- 
tization plateaux, we now present a few numerical results for correlation functions at p — 2. 

Is is technically useful to consider only objects which respect the decomposition of the 
Hilbert space according to symmetries of the Hamiltonian. We therefore define averages of 
an operator A x as 

(A) = !i (Mk) I Ei=i A X0 \Mk)) for M 0, vr, 

\^(Mk)\E L X0= AA X0 + A- X0 )\Mk)) fork = 0,n, 

where \ipo(k)) is the groundstate with momentum k. An additional advantage of this def- 
inition is that oscillations originating from the modulation of t(x) are smoothed by taken 
averages of the up to p = 2 correlation functions at a given distance. 

The connected correlation function of two quantities A and B is defined as 

C A , B (x) = (Al 0+X B X0 ) - (4 o ) (B xo ) . (4.2) 

Of particular interest are the diagonal components of the superconducting order parameter 
( p.4C| ) since quasi-long-range order is expected for one of them. 

Numerical results for correlation functions on an L = 18 system at the plateau with 
m = 1 — n are shown in Fig. f|. Characteristic oscillations are observed in the density- 
density and electron-electron correlation functions. This and the finite system size make 
a detailed analysis of the asymptotics difficult. Nevertheless, one observes that correlation 
functions containing up electrons decay faster than the corresponding ones containing only 
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FIG. 4. Correlation functions for the dimerized chain (p = 2) with U = 3t, t' = 0.7t at 
L = 18, n = 2/3 and m = 1/3. Panel a) shows density-density correlations C nain/3 (x), panel b) 
electron-electron correlations C CatCf3 {x) and panel c) superconducting correlations C^ accj ^„g{x). 
The symbols are for a = (3 =j (boxes), a = f3 =[ (filled diamonds) and a =|, =| ('x'); the lines 
are guides to the eye. 

down electrons (the latter may still be smaller in absolute value due to a smaller overall 
prefactor). In fact, all correlation functions shown in Fig. f| are very similar to those obtained 
in the non-interacting situation (U = 0) at the same L. We therefore interpret our results as 
support for exponential decay of all correlation functions containing j operators and power- 
law decay for those containing only j operators, as is expected according to the analysis of 
section |T[ In particular, these numerical results are compatible with quasi-long-range order 
in A^j on the m = 1 — n plateau. 

V. DISCUSSION AND CONCLUSIONS 

We have shown that Hubbard chains with periodic hopping or on-site energy present a 
rich structure of magnetization plateaux. More precisely, for a periodicity p, we obtain the 
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conditions (|1.2|) for the appearance of the plateaux. If both conditions are simultaneously 
satisfied, both spin and charge degrees of freedom are massive. When the combination 
pn G Z of these two conditions is satisfied, a charge gap opens irrespective of the value of 
m. Finally, if just one of the conditions ( |1.2|) is satisfied, a magnetization plateau appears 
if the total filling n is kept fixed. This result has been shown first by means of bosonization 
techniques, valid in the regime where the differences in the modulation amplitudes 5 in 
( |2.15| ) and 5/x in ( |2.20| ) are small and for arbitrary values of U. We have then shown that 
these results are confirmed by standard quantum mechanical arguments valid for small U 
and arbitrary p-merization strength and provided an expression for the width of the plateau 
to first order in U. We finally showed explicitly such plateaux in finite size systems by means 
of Lanczos diagonalization. 

The combination of a gap which can be attributed to the spin degrees of freedom and 
gapless (charge) modes prompted us to look for superconducting correlations. Indeed, we 
found quasi-long-range order in one component of the superconducting order parameter in 
those cases where only one of the conditions (|1.2|) is satisfied. 

In |H| it has been pointed out that the fully gapped situations can be most easily 
understood in the limit if — (the same argument applies also for Sa = oo). Then the chain 
effectively decomposes into clusters of p sites, whose magnetization curves are obviously 
staircase-like. The charge gap at pn G Z can also be understood in the limit of if <C t 
(5fi 3> /i), one just needs to generalize the mapping of a quarter-filled dimerized chain to an 
effective half-filled homogenous chain |3(J tom^O and commensurate filling at general p. 
Finally, this mapping can also be adapted to provide a further complementary argument for 
the existence of the doping-dependent magnetization plateaux. Again, to first order in t', 
an effective Hamiltonian can be found in the regime of strong p-merization, i.e. in the limit 
t' <^t,U (resp. 5u ^> t, U) for the case of modulation of the hopping amplitude (respectively 
of the on-site energy). When only one of the conditions (|1.2p is satisfied, which amounts to 
a condition on the filling of spin-up or spin-down bands, this effective Hamiltonian acquires 
a gap in the spin sector thus leading to a doping-dependent magnetization plateau. 

We would like to emphasize that such irrational plateaux are not present in systems 
where the doping is not fixed. Moreover, due to the remaining massless mode on such a 
plateau, the thermodynamical behavior of the system retains some particularities of a gapless 
system, such as a specific heat vanishing linearly as T — > 0. An important feature is that the 
value of the magnetization m on the plateaux at fixed n depends continuously on doping n. 
Analogous situations encourage us to believe that this scenario is generic in doped 

systems. Doping could therefore be used as a tool to study experimentally irrational plateaux 
in systems whose half-filled parent compounds exhibit plateaux only at prohibitively high 
magnetic fields. A natural candidate are ladders systems [[14]] where doping can indeed be 
controlled. Theoretical results on doping-dependent magnetization plateaux in Hubbard 
ladders will be reported elsewhere p3| . 

There are also natural problems for further study in the case of modulated chains. For 
example, the large-?/ limit of the Hubbard model leads to the t — J model. As a check of 
the generality of our results, one could therefore investigate the t — J model which at half 
filling would then be exactly the situation studied in ||. Due to the reduced Hilbert space, 
the t — J model would be particularly well suited for further numerical checks. Another 
problem to be addressed is the universality class of the transitions associated to the corners 
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of a plateau. In the case of the BA solvable model | 16|| ) it was found that the presence of a 
massless mode on a doping- dependent plateau may modify the universality class of part of 
these transitions - a fact which would be interesting to investigate also in the present model. 
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APPENDIX A: CONVENTIONS 

In this short appendix, we define our conventions and notations. The continuum fermion 
operators read 

j> a (x) = e- ik ^ R>a (x) + e ik ^ L>a (x) , (Al) 
tPl(x) = e ik ^ x ^ Ra (x) + e- ifc ^Vl,aO) . (A2) 

Using standard bosonization rules we have: 

iP R>a (x) = —^e^*-^ + . . . , (A3) 
\/27ra 

i/; L>a (x) = -^ =e ~ i ^LA^) + ... , (A4) 
V27ra 

where a is the lattice constant and (j)R t L,a are the chiral components of two real bosonic 
fields: 

4> a (x) = (f)R, a (x) + 4>L,a( X ) > ( A5 ) 

whose dual fields are defined by 

9 a (x) = <f)R ia {x) - <pL,a( X ) ■ ( A6 ) 

The dots in eqs. ( A"3) and (Al) stand for higher order terms, due to the curvature of the 



dispersion relation, which are discussed in the next appendix. The up and down Fermi 
momenta are related to filling and magnetization: 

k + = kp,] + kF,\ = nn ; k_ = kp,] — = vrm , (A7) 

where 
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n = jC^n^a) , m= j(^2s z x ) = -C£2n xA -n xd ) , (A8) 

x,a x x,a 

n x , a = c x a c x,a and L is the number of sites. Note that our definition of m (which is the one 
used for the XX Z chains in |JJ]) differs by a factor of 2 from the one of Frahm and Korepin 
23. 



APPENDIX B: FERMION FIELD OPERATOR 

In this appendix we discuss the bosonization of the fermion operator in the Hubbard 
model in a magnetic field starting from the exact BA solution. According to Frahm and 
Korepin ||23|| , the long-distance asymptotics of zero-temperature correlation functions of 



physical fields is in general a sum of terms of the form 

exp(— i2D c kp]x) exp[— i(D c + D s )kpix] 
(x — iv c r) 2A c (x + iv c r) 2Ac (x — w s t) 2A ^ (x + iv s r) 2As 

where the scaling dimensions Af s are given by: 



(Bl) 



A + / Z SS AN C - Z cs AN s y , , , 

2Af = (Z CC D C + Z SC D S ± 1\ + Nf , (B2) 

2A± = ( Z CS D C + Z SS D S ± cc 2 s detz sc c -\ + Nf . (B3) 

AN C>S , D c s , Nf s are the quantum numbers characterizing the low energy excitations. AiV c 
and AN S are integers denoting the number of electrons and down spins with respect to the 
groundstate and are fixed by the correlator under consideration. The summation runs over 
all integers or half integers D c s (depending on the parity of AA" C , AN S ) and on positive 
integers Nf, Nf. 

By analysing the leading contributions to the fermion two-point correlator, one can write 
down, after some algebra, the bosonized fermion operator: 



,0 _ Q -ik Fl x e iV^^>Ri(x) _|_ r2 Q -i2k F] x g^A-n^ _|_ ^ Q i2k F} x g 



+r 4 e~ i2kpiX e i ^*' /,i + r 5 e ~ i2k + x ^^^(h+^i) _|_ 



_j_ Q ik Fl x e -«v / 4^0L|(^)^ 1 _|_ [ 2 e i2k F fX g-iVSn^ _j_ ^ Q -i2k F ^x ^VOt^ 

+l 4 e i2kp i x e-^'t'i + l 5 e i2k + x e -iV^(H+<t>i) + . . .) (B4) 

where r iy ^ are unknown numerical constants. Notice that at U = 0, h — all these 
constants vanish except n = h — l/\ / 2na. At h = 0, the scaling dimensions of the different 
contributions in eq. ((B4|) are known from BA for arbitrary repulsion U and density n ^ 1. 
It follows that it is sufficient to retain only the following terms 



,0 — e ~ik Fi x e iV^(t>m(x)^ r ^ _|_ ^ e -i2k F ^x e iV^<t>t _|_ ^ Q i2k n x ^-i^^ -f . . .) 

_|_ e * fc F|Z e ~iV^<l>Ll(3:) n _|_ ^ Q i2k F] x g-iv^^t _|_ / 3 e ~ i2k F}X e iV^c/>t _)____) (35) 
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The expression for ijj^ can be easily obtained by exchanging j and |, with the numerical 
constants generically different. 

Using this expression for ip^ and if>^, one obtains: 



Vwj. = const d x (f>Ri + const d x (f>Li + 2rih sin(2&Fjx — \/4tt0;) 
+2(nl 2 + r 2 li) sin[2£; + :r - v^Ot + 4>i)\ 

-2(n/ 3 + r 3 h) sin[2A;_a; - Vtofa - + ■ ■ ■ (B6) 



and 

V'j^'t = c 07 ^ 9 x (pn^ + const 9^0^ + sin(2/c^x — y/Air^) 
+2(r[l' 2 + r' 2 l[) sin[2A; + a; - V^F(0 T + 0J] 

+2(r / / 3 + 7#i) sin[2A;_x - V^F(0 T " + ■ ■ ■ ( B ?) 

where we assumed the constants r, / to be real. Otherwise, the only modification would 
consist in shifts of the arguments of the sines or cosines by unknown constant phases. 

Now, assuming the constants to be equal for up and down fields, and adding eqs. 
and (|B7|) we obtain for the bosonized density operator 

p = ^j^T + = const d x ((fii + (f>i) + Ar xli sin[A; + x - Vn(4>^ + 0J] 



x cos[A;_x - v^Ot - <Pi)] + 4(nZ 2 + r 2 h) sin[2fc+a; - V47r(0 T + + • • • (B8) 

Substituting finally (fOl) , one obtains ( p.!4|) . 

Similarly, the difference of eqs. (|B~6|) and (|B~T|) yields the S z operator: 



2S Z = iptip-\ — ip\ipi = const d x (<fi<[ — 4>i) + Ar\li cos[k + x — \/vr(0t + 4>i)] 

x sin[A;_x - ^(0? _ <t>l)] ~ 4 ( r i^3 + r 3^i) sin[2A;_x - \/47r(0 T ~ + • • • ( B9 ) 



Notice the last term in the S z operator. In the usual treatments, (see for example ||26||), this 
term does not appear. As it is obvious from eq. (|B9| ), this term would be absent if we retain 
only the first two terms in eq. ( |B5[ ) or if T\ = li and r 3 = — 1 3 . 

The assumptions on the constants r i: l i} r ■, l[ to be real and equal for up and down fields 
are supported by operator product expansion computations of the original free fermion 
operator with the perturbing Umklapp operator of the Hubbard Hamiltonian, 



cos[2A; + x - V4vr(0 T + <p{)} . (BIO) 

These computations also show that, at lowest order, the constants r 2 , r 3 , / 2 , ^3 are linear in 
U. 
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